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Abstract. In this paper we provide a statistical analysis 
of the parameter-free method often used in weak lensing 
mass reconstructions. It is found that a proper assessment 
of the errors involved in such a non-local analysis requires 
the study of the relevant two-point correlation functions. 
After calculating the two-point correlation function for the 
reduced shear, we determine the expected error on the 
inferred mass distribution and on other related quanti- 
ties, such as the total mass, and derive the error power 
spectrum. This allows us to optimize the reconstruction 
method, with respect to the kernel used in the inversion 
procedure. The optimization process can also be carried 
out on the basis of a variational principle. In particular, 
we find that curl-free kernels are bound to lead to more 
accurate mass reconstructions. Our analytical results clar- 
ify the arguments and the numerical simulations by Seitz 
& Schneider (1996). 



Key words: gravitational lensing 
ies: clustering 



dark matter - galax- 



1. Introduction 

One of the most interesting applications of gravitational 
lenses is the determination of the projected mass distri- 
bution from weak lensing observations. As noted, among 
others, by Webster (1985), the mean orientation of a large 
number of distant galaxies gives a measure of the shear 
associated with the lens. The observed shear can then be 
used to derive the two-dimensional mass distribution of 
the lens responsible for the deformation induced on the 
background. This last step can be carried out in two dif- 
ferent ways. The easier route is to use a specific model 
for the lens with a number of free parameters that will 
be determined by a comparison between the observed and 
the predicted shear (see, e.g., Kneib et al. 1996). A more 
general procedure is the so called "parameter-free recon- 
struction" (Kaiser & Squires 1993; see also Bartelmann et 
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al. 1996). In this latter method the mass distribution can 
be directly determined from the shear map, provided that 
the shear is known with sufficient accuracy and detail, 
which requires the existence of a large number of source 
galaxies. 

Such reconstruction techniques are, of course, a pow- 
erful tool to study the matter distribution in clusters (see 
e.g. Tyson, Valdes, Wenk 1990, Fahlman et al. 1994, Smail 
et al. 1994) and for large scale structures. It is then im- 
portant to optimize the reconstruction process in order to 
make the best use of the observations. For this purpose, 
we have to assess the expected error of a specific recon- 
struction method, which is the main goal of the present 
paper. 

In this article we focus our attention on the parameter- 
free method, mainly because this is more general and does 
not depend on the particular lens under consideration. In 
a previous paper (Lombardi & Bertin 1998, hereafter Pa- 
per I) we have provided expressions for the error involved 
in the local measurements of the shear (or the reduced 
shear) of the lens as a function of the parameters charac- 
terizing the distribution of source galaxies. Here we extend 
the statistical analysis to the inferred global mass distri- 
bution. 

The text is organized as follows. In Sect. 2 we introduce 
the spatial weight function and we briefly describe various 
reconstruction methods used to infer the lens mass distri- 
bution. In Sect. 3 we calculate the expected error on the 
measured shear in the regime of weak lensing and in the 
more general case as a function of position in a given field 
of the sky; here the formulae of Paper I are generalized to 
the two-point correlation function for the shear map (see 
Eq. (25)). This important result is then used in Sect. 4 
to calculate the expected errors on the mass distribution 
associated with the various reconstruction methods. The 
results are then compared, in Sect. 5, to the simulations 
by Seitz & Schneider (1996). 

The main result of the paper is contained in Eq. (35) 
(together with Eq. (26)) that describes the two-point cor- 
relation function for the mass density k obtained from 
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weak lensing analysis. This proves that, in order to opti- 
mize the reconstruction process for observations in a fi- 
nite area of the sky, a curl-free kernel should be used (see 
Eq. (36)). This behavior is confirmed by numerical simu- 
lations. 



2. From the shear map to the mass distribution 

We consider a field of the sky with N source galaxies lo- 
cated at 0^ and characterized by observed quadrupole 
and ellipticity ( see Appendix A for a summary 
of the adopted notation) . Here we suppose that the galax- 
ies are observed inside a field £1 of area A, with mean 
spatial density equal to p — N/A. In the rest of the paper 
we will reserve the term "weak lensing" to the limit of 
small lens density n -C 1, i.e. 7 ~ g. 

2.1. Spatial weight function 

Source galaxies located close to a given position will bet- 
ter constrain the value of the reduced shear g{0) at such 
location. In order to describe this effect, we may thus in- 
troduce a suitable weight function W(9, 9'). The first ar- 
gument of the weight function, 0, represents the point 
of the sky under consideration and for which we want 
to measure the shear g(9), while the second argument 
0' represents the location of one observed galaxy. The 
weight function should penalize galaxies far from 9, i.e. 
W(9, + t?) should decrease for increasing ||i?||. Some ad- 
ditional "natural" conditions can be given to further char- 
acterize a specific choice of weight function and these are 
most convenient when applied (beginning with Eq. (22)) 
to a spatially continuous distribution of source galaxies 
with density p. Here we list a few possible assumptions, 
where the first is obviously the least restrictive: 

1. The weight function is even with respect to i9, i.e. 



W(6,6 + 0) = W{9,9-ti) . 



(1) 



2. The weight function is said invariant upon transla- 
tions, if it is even (see above) and if it depends only 
on the difference 9 — 9': 



W{9,9') =W{9-9') . 



(2) 



3. One natural choice is that of a Gaussian dependent 
only on the distance ||0 — 0'||, 



W{9, 9') = o _ 2 exp 



2ir pa 



if 



(3) 



where the angular scale aw should be sufficiently large 
to ensure the presence of an adequate number of galax- 
ies in a disk of radius aw centered on the generic point 
0. 

The value of the weight function at a given point 0', of 
course, has no particular meaning: only relative values are 



significant. Indeed all the following results can be shown 
to be unaffected if we merely multiply the weight function 
by a constant. Thus, we may always choose a normalized 
weight function, so that 



N 



^ W(0,0 (n) ) = 1 



(4) 

n=l 

for every 0. Of course, such a normalization will remove 
the translation invariance property, if initially present 
in W. Still, the invariance may be retained when we 
will move to a continuous description (see comment af- 
ter Eq. (22)). 

The spatial weight function operates much like the 
"shape" weight functions considered in Paper I (see 
Eqs. (23) and (21) there). In particular, using the isotropy 
condition, we can obtain the shear map g{9) either from 



N 



]T^(0,0("V( x ("), 5 (0)) = 



n=l 

or from 

N 



W yV(M (n) )Q (n) >sW =°- 



71=1 



(5) 



(6) 



In Paper I we discussed the different merits of the two 
options. In the limit of "sharp" distributions for the source 
galaxies (c <C 1), they both lead to the same determination 
of the true reduced shear go (apart from the ambiguity 
associated with the g 1— ► 1/g* invariance). 

As we will see the angular scale of the weight function 
W, i.e. the diameter of the set where W(9, 9') is signifi- 
cantly different from zero, determines a lower bound for 
the smallest details shown in the reconstructed map k(9). 
For example, in the weak lensing limit and for a weight 
function invariant upon translations, the mean value of 
the measured density k is related to the true density kq 
through the expression (see Appendix C) 



(k(0)) = p J W{9 - 9')k q (9') dV . 



(7) 



This relation, similar to that found when an image is de- 
graded by a PSF, suggests the possible application of de- 
convolution techniques (see Lucy 1994) to the present con- 
text. 

2.2. Weak lensing regime 

We first consider the case where f2 is very large, so that 
the field is identified with the whole plane (the effects of 
the boundaries will be discussed soon, in Sect. 2.4). 

There are basically two ways to reconstruct the mass 
distribution n{9) from the shear map g{9) (Seitz & Schnei- 
der 1996). The first, more natural method is based on the 
integral relation 



k(0) = J A(0-0'H(0')cfV 



(8) 
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(9) 



where the kernel Di is given by (Kaiser & Squires 1993) 

D 1 (6)\ J_ 
D 2 (9) ) 7r||0|| 4 

In the weak lensing limit 7 ~ g, and thus the reduced 
shear map can be used directly in Eq. (8) to derive n(9). 
Note that the inverse relation holds with the same kernel 



9\ -el 
26182 



7 ,(0) = J D l {9-9')n{9')A 2 d l ■ 



(10) 



This relation will turn out to be useful in Appendix C. 

A second possibility, which can be proved to be math- 
ematically equivalent to the first, is based on the exact 
relation (see Kaiser 1995) 



Vk = u(6) 



71,1 + 72,2 
72,1 - 71,2 



(11) 



which is a direct consequence of the thin lens equations. 
Here 7^ = dji/dOj. By analogy with the condition used 
to derive Eq. (8), if we assume that n(9) vanishes for large 
values of ||0||, Eq. (11) can be inverted to give 



n{9) = j Hf s {9 - 9') Ui {9') dV , 



with the kernel 



(12) 



(13) 



In Eq. (12) the shear map enters through the vector u, 
which, in the weak lensing limit, involves the derivatives 
of g{9). This second method thus introduces undesired 
differentiations, but it has the advantage that it is more 
easily generalized to include the effects of the boundaries 
(see Sect. 2.4 below). 

2.3. The general case 

When the lens is not weak, Eq. (8) 

k{9) = j D t {9-9'){l-K{9'))g t (9')d 2 8' (14) 

can be solved by iteration (Seitz & Schneider 1995). 

The second method, related to Eq. (12), has been gen- 
eralized by Kaiser (1995) for strong lenses. If we introduce 
k(9) = ln(l — k(9)) and the new vector 



\g\ 2 



1 + .91 

92 



92 

1-51 



.91,1 + #2,2 
.92,1 - ffl,2 



(15) 



then it is possible to show that the relation Vk(0) = u(9) 
holds. As a result k can be obtained from u via the same 
integral equation (12) used earlier. The fact that k is de- 
termined only up to a constant here translates into a non- 
trivial invariance for the density distribution k(9), under 
the transformation (see Schneider & Seitz 1995) 

k(0) 1— > (1 - C)k(0) + C , (16) 

consistent with Eq. (14). 



2.4- Effect of the boundaries 

The methods described so far assume an infinite domain of 
integration. In practice, one can measure the shear only in 
a finite area (e.g. the CCD area), which is often small com- 
pared to the angular size of the lensing cluster. Therefore, 
the relations given earlier should be properly modified. 

We briefly noted that the second method is better 
suited for the purpose. In the following, for simplicity, we 
consider only the weak lensing limit, but the equations 
that we will provide can be easily extended to the general 
case by replacing (it, k) \—> (u, k). The relations suggested 
by Seitz & Schneider (1996) for mass reconstruction in a 
field £1 of finite area A are of the form 



k(6)-k = [ G l {9,9')u l {9')& 2 6' 
Jn 



(17) 



Here R is a constant representing the average of k, while 
G is a suitable kernel. The kernel is chosen so as to give 
the correct mass distribution if u could be measured with 
no errors (see Eq. (C14)). There is however some freedom 
left in the choice of the kernel, mainly because it returns 
a scalar field (k) from a vector field (it). This freedom will 
be further discussed later on. One interesting kernel, called 
noise filtering, has been introduced by Seitz & Schneider 
(1996) 

H ss (8,9') = -V e ^ ss (0,0') , (18) 

where H ss is the solution of Neumann's boundary problem 
(n is the unit vector orthogonal to dCl) 



V 2 g ,H ss (9,9') =5(9-9') - 
dH ss (9,9' 



89' 



-rn{9') = 



A ' 

W e on . 



(19) 
(20) 



The term related to the area A ensures the proper appli- 
cability of the Gauss theorem. Note that the kernel H ss 
has vanishing curl, i.e. 



V e < AH bb {9,9') =0 . 



(21) 



3. Measurements of the reduced shear map and of 
the two-point correlation function 

In this section we will give an expression for the reduced 
shear map measured using Eqs. (5) or (6). In Paper I we 
have calculated the statistics associated with a local shear 
measurement under the hypothesis that the probability 
distribution for the source ellipticity x s is sharp, i.e. most 
source galaxies are nearly round. Now we consider situa- 
tions where the reduced shear is a function of the position 
0, but we assume that g(9), a smooth function of 0, does 
not change significantly on the angular scale aw of the 
weight function W(9,9'). An important new aspect of the 
analysis that has to be addressed here, in view of the goal 
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of determining the error on the reconstructed mass, is the 
calculation of the two-point correlation function for the 
shear map. 

So far we have considered source galaxies with ran- 
dom orientation but with fixed position on the sky (corre- 
sponding to 0^ on the observer's sky). It is interesting to 
average all results by assuming that galaxies have random 
positions. The result of the average can be approximated 
by considering a continuous distribution of galaxies with 
density p (number of galaxies per steradian). This leads 
us to ignore, for the moment, the effects of Poisson noise 
associated with the finite number of source galaxies (fur- 
ther comments are given at the end of Sect. 4.1). Following 
Seitz & Schneider (1996), we consider a homogeneous dis- 
tribution of galaxies in 0, i.e. in the observer's plane. If 
p is independent of 0, we may change summations with 
integrals using the rule J2 n p J d 2 8' . Here, to simplify 
the derivations and the discussion, we adopt, for every 0, 
the normalization 



p J w(e,e')d 2 e' = 



i 



(22) 



for the weight function. For an infinite field Q, or for the 
case described in Fig. 1, this normalization does not break 
the translation invariance of W. Then as shown in Ap- 
pendix B, the relation between expected and true value of 
g, corresponding to Eq. (7), is 



(9(0)) 



w(e,e')g (e')d 2 9' . 



(23) 



As is intuitive, "near" galaxies give the most important 
contribution to the measured value of g. 

The correct generalization of the covariance matrix 
Covij(g) when g is a function of the position is a two- 
point correlation function: 



CovyfoM') 



((g i (e)-^m)(g j (e , )- 9i {e'))). 



(24) 



In the weak lensing limit Eq. (25) then reduces to 
Cov ij ( 7 ;0,0') = ^ / W(0,0")W(0',0")d 2 8" 



4 

cS 



w 



\0 - 



'112 



(26) 



Note that the knowledge of the "diagonal" values Covy(g; 0, 
is not sufficient to calculate the error on other variables, 
such as the density distribution k, determined from g (cf. 
Eq. (32) and Eq. (35)). 

If we assume that the weight function is even (property 
1 of Sect. 2.1), then the two-point correlation function of 
g can be written in the simple form (sec Appendix B) 



CovyfoM') = 4(1 - I (.9(0)) | 2 ) (1 - |<.9(0')>| 2 )% x 



p J W(0,0")W(0',0")d 2 0" . 



(25) 



Here c is the covariance of the ellipticity distribution of the 
source galaxies. In this equation, as noted for Eq. (23), we 
suppose the weight function W(0, 0') to be normalized. 



The last relation holds for a Gaussian weight function of 
the form given in Eq. (3); here Cov^-^; 0, 0') is a sim- 
ple Gaussian with variance 2er^ and depends only on 

||0 - 6'\\. The variance ((7,(0) - 7i(0)) 2 ) of 7,(0) is sim- 
ply Covj,( 7 ;0, 0) = c/(16npa w ) (without summation on 
i) and thus increases if aw decreases. This behavior can be 
explained by considering that the number of galaxies used 
for a single point is of the order of pa w . Notice also that 
aw sets the scale length of the covariance of 7: measure- 
ments of 7(0) and 7(0') are uncorrelated if ||0— 0'|| 3> aw- 



4. Measurements of the mass distribution 

It is not difficult, at least in principle, to calculate the 
error on k(0) from the two-point correlation function of 
g. The error on k, of course, depends on the reconstruc- 
tion method used. For this reason, following Sect. 2, we 
consider different methods separately. For simplicity, we 
suppose that the weight function W is invariant upon 
translations. Moreover, we suppose that the angular scale 
of the weight function W is much smaller than the an- 
gular scale of k (i.e. the scale where n varies signifi- 
cantly). In general, if we ignore edge effects, the rela- 
tion between the error on g and that on k is given by 
Covjj(7; 0, 0') = Sij Cov(k; 0, 0'). Here we show only the 
results obtained, referring to Appendix C for a derivation. 
As in Sect. 2, we first refer to the case where is identified 
with the whole plane (finite field effects will be addressed 
in Sect. 4.3 below). 

4-1. Weak lensing regime 

In this case we can use either Eq. (8) or Eq. (12). A rather 
0)jprising result is that both methods lead to the same 
mean values and errors for k. The result for the mean 
value has already been stated in Eq. (7), i.e. the measured 
mass distribution (k(0)) is the convolution of the weight 
function W with the true mass distribution kq(0). 

For an "isolated lens" (a case where n is taken to vanish 
outside a certain domain Qi n ) the ambiguity associated 
with Eq. (16) is resolved and the concept of total mass of 
the lens becomes meaningful. In the weak lensing limit, 
from any reconstructed k(0) one can in principle accept 
also Kc(0) = «(0) + C . Now if we know that the density 
vanishes outside Sl\ n , the constant C can be determined 
by requiring 



/ 



d 2 0K C {0) = 



(27) 



M. Lombardi & G. Bertin: Accurate mass from weak lensing 



5 



where f2 out = fi \ Oj n is the part of the field not contained 
in fi in . Therefore, the appropriate density to be used is 



kc(0) = k{0) 



1 



A>ut Jn c 



n(e')d 2 e' 



(28) 



The associated total mass is 
M = [ Kc (8)d 2 9= f K C (9)d 2 e 



/Q, n 

Therefore: 



K (e)d 2 e 



(M)= / (n{0))d 2 



Hi, 

A- u 

4 Q ut jo, 



A\r 



= / K (0)d 2 6 = M 

Jn ia 



n{0)d 2 9 . 



( K (e))d 2 



(29) 



(30) 



where M is the true mass of the lens. In other words, the 
smoothing effect associated with W does not change the 
measured total mass M of the lens. 

The covariance of the lens distribution k can be shown 
to be equal to (for both Eqs. (5) and (6)): 



we expect J\f = pL 2 /s 2 galaxies per patch. In this case the 
expected variance of 7 is (see Paper I) c/4Af, and the vari- 
ance of k is of the same order of magnitude. The expected 
variance of M is then (c/4/V) (L 2 /s 2 ) s 2 , where the first 
factor is the variance of n in every patch, the second factor 
is the area of every patch (the square is necessary because 
we are dealing with variances), and the third factor arises 
because we must add s 2 independent variables. The final 
result for the variance of M is cL 2 /4p, exactly as stated 
byEq. (32). 

Here we may come back to the issue of the Poisson 
noise, only mentioned at the beginning of Sect. 3. Strictly 
speaking, the relation Af = pL 2 / s 2 in the previous para- 
graph should be replaced by (AO = pL 2 /s 2 , with Af follow- 
ing a Poisson distribution. We now consider the variance 
of 7 as a function of Af. An estimate of the effect of the 
Poisson noise can be obtained in the limit (M) !3> 1 by 
expanding 



Var( 7 ,A0 



4<A0 



Af - (AT) (Af - (AT) 



(AO 



(AO 



Averaging over the ensemble thus yields 



Cov(k; 0, Q')= C -j J W{0, 0")W{0', 0") d 2 6" . (31) (Var( 7 , A/")) 



4<A0 



1 



1 



•(33) 



(34) 



In comparing Eq. (31) to Eq. (26), one should note that 
the similarity of results refers statistically to average er- 
rors, but not to the individual errors of one reconstruction. 
The variance in the measure of the total mass is the inte- 
gral of the covariance of n: 

Var(M) = f d 2 6 [ d 2 6' Cov(k; 0, 0') 

I d 2 e f d 2 e' Cov( K ; 0, 0') 



A\ 



lout 



cA 



(32) 



where, we recall, A is the area used.* Obviously, the latter 
approximate expression holds when A; n /^4out <C 1. Curi- 
ously, this result does not depend explicitly on the weight 
function W . The derivation given in Appendix C assumes 
that the weight function is of the form of Eq. (3), but a 
similar expression for the variance of M is expected to 
hold in the more general case. 

The results of this subsection can be clarified by a sim- 
ple example. Instead of introducing the weight function 
W, we consider the unweighted Eqs. (5) or (6) on small 
patches of the sky. For simplicity, we refer to a square set 
ft of length L divided into s 2 equal square patches: thus 



* In principle, all integrations should be performed in 
the whole plane. However, here we suppose to perform 
integrations over a domain 17 of area A large with respect 



to a 



w ■ 



The effect of the Poisson noise is here contained in the sec- 
ond term in brackets, which is negligibly small. Additional 
discussion is postponed to the end of Appendix B. 

4-2. The general case 

The situation is, in principle, quite similar to the weak 
lensing limit, but, in practice, the calculations are much 
more difficult. If the angular scale of k is much greater 
than the angular scale of W, then we can prove that the 
mean value of the measured mass distribution given by 
Eq. (7) holds unchanged. 

Difficulties in the calculation of the covariance of k 
mainly arise from the form of the covariance of g given 
by Eq. (25), because of the dependence on and 0' of 
the first factor. However, if the lens has |<7o(#)| < 1, the 
covariance given by Eq. (25) is smaller than that of the 
weak lensing limit of Eq. (26) and thus we can consider 
all the results given in the weak lensing regime as upper 
limits for the errors. 

4-3. Edge effects 

Finite boundaries introduce interesting effects, and make 
the errors depend on the kernel G used in Eq. (17). For 
simplicity we take two different sets for (see Fig. 1). The 
first set is fi', i.e. the observation area that includes all the 
lensed galaxies used in the reconstruction. The second set 
is fi C ft', i.e. the set where we measure g{0)- We suppose 
that every point in f2 has a neighborhood with radius of 
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Fig. 1. Sketch of the observation area used in the mass recon- 
struction. 

the order of the angular scale of W completely enclosed 
in Q' . This assumption greatly simplifies calculations and 
does not have major practical consequences, except that it 
leads to discarding a small strip O' \ O around the bound- 
ary <90' of 0'. With this hypothesis the expected mea- 
sured mass distribution is again given by Eq. (7) as long 
as <G O. However this is strictly true only if we choose 
correctly the mean mass distribution R (see Eq. (17)). 

In general, the covariance depends on the kernel G 
used, and in particular on its divergence-free component 
(see Appendix C). [We recall that a vector field G(6,9') 
can be decomposed as G = G' + G", where G' has van- 
ishing curl and G" has vanishing divergence (as usual, in 
the above notation and decomposition, the emphasis is on 
the variable 9', since 9 is taken to be fixed).] The result is 

Cov(k;0,0) = Cov( 7 ; 0, </>) 

+ / d 2 9'[ d 2 0'G;'(0,0')Cov( M ;0',0')G"(0,0') , (35) 
Jn Jn 

where u is the quantity defined in Eq. (11). The first term 
is clearly independent of the kernel G used, while the sec- 
ond term can be shown to be positive definite, i.e. 



Cov(k; 9,9) > Cov( 7 ;0,0) . 



(36) 



Thus the error on n is minimized if a curl-free kernel G 
is used. This suggests that only curl-free kernels should be 
used in weak lensing reconstructions. In fact, the kernels 
so far judged to be "good" by means of simulations, all 
have vanishing curl (see Sect. 5). For a curl-free kernel, 
such as the noise filtering kernel given in Eqs. (18-20), 
the result is independent of the kernel used and of the set 
0. 

We now investigate the class of kernels G(9,9') that 
satisfy the following properties: 



i. G inverts Eq. (11) when u{9) is measured with no 
error; 

ii. G is curl-free. 



From the second property we can write 
G{9,9') = -V e ,g(9,9') . 
Thus, if u{9') = Vk(0') wc find 

k{9) -r = - f Vg{9, 9') ■ Vk(0') dV 
Jn 

K (9')V 2 g(9,9 , )d 2 0' 



(37) 



/ ■ 

Jon 



K (9')Vg(9,9') -nd6' . 



(38) 



As in Sect. 2, n is the unit vector orthogonal to <90. The 
last relation shows that, in order to satisfy point i, we 
must have 



V 2 £(0, 9') = 5(9 - 9') - 



(39) 



If k(9) is not known on the boundary of O, we should also 
consider: 



dg{9,9') 



Hi =0 



V0' G dfl 



(40) 



In this case the kernel to be used is simply the one ob- 
tained from Q = H ss , i.e. G = H ss (cfr. Eqs. (18-21)). 
Otherwise, Eq. (40) is to be dropped and the kernel g is 
determined up to a term C, where £ is a harmonic func- 
tion (V 2 £ = 0). This free function can be used to dispose 
of the contribution that would arise from the boundary 
term in Eq. (38). 

As, in general, the measured u field is not curl-free, 
the inversion can only be approximate. The best inversion 
can thus be found by searching for the function n(0) that 
minimizes the functional 



S 



|Vk(0) -u{9)\\ d 2 e . 



(41) 



If we vary the distribution n k + 5k, the functional 
would in general change to S + SS, with 



5S 



= 2 f V(5k) ■ (VK-u)d 2 6 
Jn 

= -2 f 5k(\7 2 k-X7 ■u)d 2 9 

+ 2 5k(Vk-u) - ndO . 
Jan 



(42) 



By setting 5S — we readily find the associated Euler- 
Lagrange equation 



V 2 k(6») = V • u{9) . 

This equation should be supplemented by 
Vk(0) • n = u{9) n V0 e dQ , 



(43) 
(44) 
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unless we fix the value of k on dil, so that the boundary 
term in Eq. (42) vanishes because 8k = on dft. Equa- 
tions (43) and (44) define a Neumann boundary problem 
equivalent to Eqs. (39) and (40) above, in the sense that 
G is precisely the Green function associated with it. This 
clarifies the interesting properties of H ss . 

4-4- Power spectrum 

In order to express in a simple manner the errors involved 
in the reconstruction process, Seitz & Schneider (1996) 
introduce a "power spectrum" P(k). In the weak lensing 
limit, their definitions are 

A(fc) = ±J e k e ( K (6) - ko(0)) d 2 6 , (45) 

P(*) = (A(fc)A*(fe)), (46) 

where the mean in Eq. (46) is also over the various direc- 
tions of k. As a result, the power spectrum P(k) is simply 
the variance of the complex map A(fc), i.e. the Fourier 
transform of the reconstruction error. Thus, for example, 
the value of P(0) is proportional to Var(M), the error on 
the total mass, while its behavior for larger values of k is 
related to the angular scale of the weight function used. 

Within our framework it is not difficult to evaluate 
the relevant power spectrum. A simple calculation (see 
Appendix D) for a Gaussian weight function gives 

P( fe )~_£_exp[-^A ! 2 ] , (47) 

i.e. a simple Gaussian with variance l/2a^. One might 
anticipate a significant contribution to the power spec- 
trum coming from the error associated with the difference 
between («) and kq, but we argue in Appendix D that 
such contribution is negligible in the weak lensing limit. 

5. Comparison with numerical simulations 

In this section we compare our predictions with the re- 
sults obtained by Seitz & Schneider (1996) from numeri- 
cal simulations. Simulations start by defining a lens mass 
distribution and a random sample of source galaxies. Each 
galaxy is traced to the lens plane and the reduced shear 
g is then calculated from the observed ellipticities using 
Eq. (5). Finally the shear map is inverted into the lens 
mass distribution n using various methods. For Seitz & 
Schneider take a square 7.5' x 7.5'. Source galaxies have 
random orientations and their ellipticities follow truncated 
Gaussian distributions. Simulations have been performed 
with three different variances for x s : c\ = 0.069109, 
c 2 = 0.13323 and c 3 = 0.19689. 

The reconstruction method used by Seitz & Schneider 
is similar to the one described in Sect. 2, with the following 
differences: 



i. Their weight function is not invariant upon transla- 
tions because it is a Gaussian of argument \\6 — 8'\\ 
with the variance depending on 6. 

ii. An outer smoothing is added to the final lens mass 
distribution. 

The first point is a device introduced in order to have 
better resolution in the stronger parts of the lens. The 
second point is used in order to have a smooth lens distri- 
bution from a discrete map of k. 

Our result (47) for the power spectrum can be easily 
generalized in order to take into account the outer smooth- 
ing: 

F(i)^exp[-(^ + ,^ 2 ] . (48) 

Here o\ is the variance associated with the outer smooth- 
ing. Note however that the expression given above does 
not take into account the variable-scale smoothing used 
by Seitz & Schneider. Even if this result has been derived 
with some approximations (weak lensing limit, large area 
A of f2, fixed inner smoothing aw), a comparison with the 
simulations shows that Eq. (48) can reproduce the main 
features of the simulated power spectrum. 

Figure 2 shows the results of simulations together with 
the power spectrum predicted by Eq. (48). Thus Eq. (48) 
underestimates the error. This difference can be attributed 
to the following factors: 

1. The weight function W considered is not precisely of 
the form of Eq. (3), because of the change of normal- 
ization near <9f2. This last factor should increase the 
variance of k near the boundary of fi (the variance of 
7, with direct influence on n, should double near a side 
of £1 and quadruple near a corner; cf. top-right frame 
of Fig. 10 in Seitz & Schneider 1996). 

2. The weight function W is not of the form of Eq. (3) 
also because of differential smoothing. 

3. The constant term R has not been estimated prop- 
erly (see first paragraph of Sect. 4.3). In principle, this 
should be traced to a counterpart in P(0), but for finite 
sets il there is an additional term in P(k). 

4. The set il is not the whole plane. 

5. The lens is not weak (see Eq. (25) and the extra con- 
tribution in Eq. (D4)). 

6. Poisson noise is associated with the finite number of 
source galaxies. 

7. The population of source galaxies is characterized by 
sizable c. 

In spite of these limitations, the general behavior of 
P(k) is reasonably well reproduced by Eq. (48). In par- 
ticular, the maximum of the simulated points corresponds 
exactly to the maximum of our theoretical curve. 
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simulations. In fact, as stated in Eq. (21), this kernel is 
curl-free. 

Another kernel considered is the "S-inversion" (see 
Schneider 1995). Simulations show that errors of the S- 
inversion are nearly the same as for the SS-inversion. The 
S-inversion operates by averaging over radial paths made 
of two segments. Inside the segments the kernel is easily 
shown to be curl-free. 

The last kernel is the so called "B-inversion" (Bartel- 
mann 1995; see also Squires & Kaiser 1996). From the 
results of the simulations it is clear that the B-invcrsion 
leads to larger errors on the map distribution. This be- 
havior is once again explained by Eq. (35), since the B- 
inversion kernel is not curl-free. Notice that in this case it 
is difficult to estimate analytically the exact error on k. 

Acknowledgements. We would like to thank Peter Schnei- 
der for several comments and suggestions that have helped 
us improve the paper. This work has been partially sup- 
ported by MURST and by ASI of Italy. 



Fig. 2. Comparison of predicted power spectrum (dashed lines) 
with measured power spectrum (solid lines) for the simulations 
made by Seitz & Schneider (1996). All frames refer to a source 
population characterized by c = 0.069109. Top frame: pA = 80, 
a w = 0.212', o s = 0.0778'. Middle frame: pA = 50, a w = 
0.177', a s = 0'. Bottom frame: pA = 50, a w = 0.240', a s = 
0.0778'. 

5.1. Curl- free kernels 

In order to check the result of Eq. (35), we have considered 
different kernels used by various authors and we have com- 
pared our predictions with other aspects of the numerical 
simulations performed by Seitz & Schneider (1996). 

The first kernel considered is the noise-filtered "SS- 
inversion" (Seitz & Schneider 1996) described above in 
Eqs. (18-20). This method has been especially designed 
to reduce the statistical errors and performs very well in 



Appendix A. Notation 

We collect here the main symbols used in this paper. 

X s Superscript s identifies source (unlensed) quanti- 
ties. 

Xi Subscripts refer to the complex representation, 
e -g- X = Xi +1X2 or, when indicated, to the vector 
representation, e.g. 9 — {61,62). 

9 s , Unlensed and observed position of a point source. 

E c Critical density. 

£(0) Projected mass distribution of the lens. 

k{0) Dimcnsionless mass distribution: k{0) = E(0/E c . 

Q s , Q Unlensed and observed quadrupole moment of an 

extended image. 
X s , X Unlensed and observed (complex) ellipticity of an 

extended image: x s = (<2u -Qh+^Qh) / (Q11 + 

Q| 2 ), and similarly for x- 
7 Complex shear. 

g Complex reduced shear: g = 7/(1 — k). 

c Covariance of the source ellipticity for an isotropic 

distribution: (XiXj) = c %- 
Id, Sij Identity matrix. 

Appendix B. The shear 

In this Appendix we derive Eq. (23) and Eq. (25) of Sect. 3. 
We assume that the reduced shear g has been measured 
through Eq. (5). Calculations based on Eq. (6) are very 
similar. As explained in Appendix A of Paper I, the mean 
value of g obeys the relation 

JV 

]T w{0, 9^) x s (<x (n) ), (g{0))) = . (Bi) 

n=l 
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Here (x^) is the mean value of \ m 0^ an d thus de- 
pends on g (9^). We now assume that go(0) does not 
change significantly on the angular scale of W(9, 9'). This 
implies that we can expand the previous equation to first 
order in g . We choose, as starting point, the value 50 = g* 
given by 



N 



71=1 

Then we find easily 



(B2) 



N 

Y, w * s (x.,(g{0))) 

71=1 



JV 



71=1 



x=x. 



dg* 



(a>(* (B) 



, (B3) 



where x* is the expected value of the ellipticity when 
the reduced shear is equal to g*. Equation (B3) has the 
obvious solution (g(0)) = g*(6), because the first term, 
X s (x* > <?*) 7 vanishes by definition and the second vanishes 
due to the choice of <?* (notice that the partial derivatives 
in the latter do not depend on 6^). When moving to a 
continuous description, we have to calculate the average 
expected value of g over all possible positions {0^} of 
the source galaxies 



(B4) 



9 {0) -JV l ' 1 "" ' [ ' flj 2 I d ^ , V) > 

N 

xY J w{e,e^)g Q {e^) 



71=1 



The weight function W(0, 6 {n) ) depends on all { (1) , (2) , . 
because of the adopted normalization. However, in the 
limit N » 1, the weight function W(Q,Q {n) ) can be con- 
sidered to depend only on 0^ alone, so that the above 
integral can be approximated by 



1 



JV 

(9(0)) - j £ 

71=1 



W(6,0<ri)go(oW 



P 



d 2 e'W(0,0')g o (0') . 



(B5) 



This proves Eq. (23). Our result simply states that the use 
of the first order expansion in g reduces every mean to a 
weighted arithmetic mean. 

Calculations for the covariance of g are much more dif- 
ficult but basically repeat those given for the unweighted 
situation in App. A.l of Paper I. In particular, if we call 
F({8 {n \x {n) },g(0)) the function defined in the l.h.s. of 
Eq. (5), we have 

Cov(g;0,0')=B- 1 (0) 



JV 



E^ (ll) (0)Cov( X ( "V ( " )T (0') 



(B-\e')y 



(B6) 



where 
B(9) = 
A< n >(0) 



dF 
dgJO) 1 
dF 



(B7) 
(B8) 



All functions have to be calculated in the mean value of 
their arguments. Some calculations then lead to 



B{6) 



1 



\(g{0))\ 



Id . 



(B9) 



The term in brackets in Eq. (B6) can be written in the 
form 



Y,A^(8)Cov(x {n) )A^ T (0') 

71=1 

JV 

= cY J W(e,e {n) )W(0',e {n) )+\meBX terms . (BIO) 

71=1 

Here "linear terms" means additional terms linear with 
respect to the quantity [go(0^) —g*\ , based on the same 
expansion defined by (B2). 

By averaging over the source positions and by moving 
to a continuous description (basically following the same 
steps indicated in (B4) and (B5)), we thus obtain Eq. (25). 
Notice that the "linear terms" in Eq. (BIO) do not give 
any contribution when averaged over the source positions. 
We stress that the results stated here are valid only if the 
weight function is even (property 1 of Sect. 2.1). 

^•paily, we point out that the approximation that takes 
'lis frorri p54) to (B5) is precisely associated with neglecting 
the Poisson noise. 



Appendix C. The lens mass distribution 

In this Appendix we will derive Eq. (7) and the results 
stated in Sect. 4, assuming a weight function invariant 
upon translations (case 2 of Sect. 2.1). 



C.2. Weak lensing limit 

Calculations in the weak lensing limit are not difficult. 
As explained in Sect. 2.2, we can use either Eq. (8) or 
Eqs. (11) and (12) to convert the reduced shear into the 
mass distribution. 

In the case of Eq. (8) we can write 



(k) = Di* (7j) = A* [pW * 7o») 



(CI) 



Here the star denotes convolution, while 7oi is the com- 
ponent i of the true shear map 70. The second step is 
justified by Eq. (23) applied in the weak lensing limit. By 
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using the associative and commutative properties of the 
convolution and by noting that k = A*7oi) we can write 

(k) = Di-kpW-k 7oi 

= pW*A*7to=pW'*Ko- (C2) 

This proves Eq. (7). 

About the covariance of n we can write 

Cov(k;6>,0)= / dV / dV Di(d - 0') x 



(C5) 



xCov( 7 ;0'-0')A(0 -</>')• (C3) 

As A and Cov(7) are even functions, this is simply a 
double convolution, and thus the result depends only on 
the difference between and <fi. Therefore we can write 
for Cov(k; 0, </>) = Cov(k; — (f>) the expression 

Cov(k) = A * Cov( 7 ) * A = Cov( 7 ) , (C4) 

i.e. Eq. (31). Here we have used again the commutative 
property of convolutions and the relation A * Di = 5 
given by Eqs. (8) and (10). [Hereafter 5 means the Dirac 
delta distribution.] 

The results are the same if we use Eqs. (11) and (12). 
In fact we can write Eq. (11) as a convolution between 7i 
and the operator 

5,i(0) 6,2(6)' 
-6, 2 (0) 5,1(0), 

where 6,i(0) = 85(0) /89i. Thus we are allowed to use 
the properties of convolutions. It is obvious then that the 
convolution with the weight function W in ( 7 ) can be 
moved to the true lens distribution kq, and we find again 
the result of Eq. (C2). The covariance matrix of u can be 
calculated using the operator (C5). As a result, we find 

Cov(u) = -V 2 Cov( 7 ) . (C6) 

We then have 

Cov(k) = -Hf s * V 2 Cov( 7 ) * Hf s 

= -V 2 (#P*ffP)*Cov( 7 ). (C7) 

A simple calculation shows that -V 2 (i? I KS * Hf s )(0) = 
V 2 (ln ||0||/27r)= 5(0). This leads again to Eq. (31). 

Now let us prove Eq. (32). From Eq. (C4) and Eq. (26) 
we find 



Cov(k) = ^-W -k W 
4 



(C8) 



In the limit A- m /A out <C 1 explained in Sect. 4.1, the main 
contribution to the variance of M derives from a double 
integration of Cov(k). A simple change of variables gives 



Var(M) = ^ J d20 J d 2 6' (W*W){0') 
cpA 



rA 

4 : W(0)W(0) = - 



(C9) 



Here W is the Fourier transform of W and the last equal- 
ity holds because of the normalization (22) of the weight 
function. 



C.3. The general case 

In the general case we restrict ourselves to estimating the 
mean value of the lens distribution because calculations for 
the covariance are too difficult. Under the hypothesis that 
the angular scale of W is much smaller than the angular 
scale of k (or g), the situation is much like that of the 
weak lensing limit. As shown in Appendix B, this basically 
implies that all averages are weighted arithmetic averages. 
Simple calculations show that we have (u) = iiQ-kpW, and 
hence (k) = kq * pW . As usual the assumed ordering of 
scale lengths leads again to Eq. (C2). 

C4- Edge effects 

For simplicity we refer to R = 0. We rewrite Eqs (17), (11) 
and (10) with a different notation 

7TQK = GiUi , (CIO) 

Ui = 3ij*7j. (CH) 
7i = A * « • (C12) 

Here ttq is the characteristic operator for the set f2: 

M)(0) = { f n {0) if f f e "' (C13) 
i. otherwise . 

Equation (CIO) is equivalent to Eq. (17) with R = if 
we redefine the kernel G(0, 0') for every and 0' so that 
G{0, 0') = if either fl or 0' £ Q. With this simple 
definition we can extend the integration domain (usually 
fl) to the whole plane. Notice that while and Di are 
used in convolutions, is a generic linear operator. From 
these equations we have 



7tqk = dui = G t (Tij ★ Dj ★ k) , 
and thus we find the identity 
GiTij ★ Dj = 7Tsi . 

Equation (12) with the new notation is 
k = H^ s * Ul . 



(C14) 



(C15) 



(C16) 



This, together with Eqs. (Cll) and (C12), gives us another 
identity: 



★ Dj -k Hf s ★ u k . 



(C17) 



Using Eqs. (CIO), (C12), and the relation (m) = pW-ku i, 
we can easily obtain the mean value for measures of k: 



(k) = Gi(ui) = GipW-kUoi 

= dTij ★ Dj k Hf s * u ok * P W 



TTnHf s k u k * pW = TT n pW * k 



(C18) 



As usual, subscript indicates the true value of a quantity. 
This equation, rewritten in the more standard notation, 
is Eq. (7) for G 0. 
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Lot us calculate the covariance of re. First of all note 
that, while Eq. (C14) implies Eq. (C15), from Eq. (C17) 
we cannot deduce that * Dj * H^ s is the identity. 
This happens because the two components of u,- L arc not 
functionally independent, as one can see from the relation 
V A u — 0. In fact, using Fourier transforms it is easy to 
prove that the operator R ik = * Dj * H^ s selects the 
curl-free component of a vector field. Its Fourier transform 
is 



Rij(k) 



k^kj 



From (CIO) we have 
Cov(re) = GjCov(u)Gj . 



(C19) 



(C20) 



Every vector field can be written as the sum of two vector 
fields, of which one is curl-free and the other is divergence- 
free. Hence, if we consider G(0, 0') a vector field with 
respect to 0' , we can write 



G(0, 6') = G'(0, 0') + G"{0, 0') 
where 

V e > AG' (9, 6') = 0, 

V > -G"(e,e') = 0. 



(C21) 

(C22) 
(C23) 



Thus G' and G" can be written as the gradient and the 
"curl" of two scalar fields: 



G'{e,e') = v e > s '{e 1 o') , 

G"{0,0') = V«,< As"{0,9') 



4(0,0') 

-s'[(0,0') 



(C24) 
(C25) 



There is some freedom in the choice of G' and G" (or 
equivalently of s' and s"). However, it is always possible 
to choose G' and G" so that they vanish for — > oo. 
With the decomposition (C21) we have GiRij = G'j and 
thus 



Cov(re) = {G t R l3 + G")Cov(u)(G k R kl + G") 



(C26) 



Recalling now the definition of Rij and using Eq. (C15) 
we find 

Cov(re) - (wnHf s + G;')Cov( u )(7r^i/f s + G") (C27) 

If the kernel G has vanishing curl, then G" = and we 
find the final result 

Cov(re; 0, 0) = J d 2 0' J dV Hf s (6 - 0') x 

x Cov(u; 9' - </>')fff S (0 - 0') for 0, e . (C28) 

In general however we must evaluate three additional 
terms. Two of them are of the form 

G';Cov{u)G' 3 = J dV J d 2 ^'G"(0,0') x 

xCov(u;6>'-0')-V^s'(0,0'). (C29) 



where we have used Eq. (C24). By the change of variable 
4>' 6' — $ and after integrating by parts we find 

G;'Cov(ii)G; = j d 2 9' J d 2 <p'V e > ■ G"{6,&) x 

x Cov(u; 0' - cf>')s"((j) - 0') = , (C30) 

where the last relation holds in virtue of Eq. (C23). We 
finally rewrite Eq. (C27) in a simplified form: 



Cov(k) = TT n H^Cov(u)TrnH^ + G;'Cov( w )G;' (C31) 

Here the first term is independent of the specific kernel 
G used, while the second term depends only on G" . As, 
by definition, Cov(m) is positive definite, the last term in 
Eq. (C31) is also positive definite. In other words if G" ^ 
the error on re will increase. 



Appendix D. Power spectrum 

The power spectrum reported in Eq. (47) can be deduced 
from the expression of the mean and covariance of the 
measured mass distribution. In particular we have 



(A,(fe)) = l^( c ^((re(0)) -n o {0))d 2 9 
Covij (A; k, k') = -L J d 2 6 £ d 2 9' (c ife e ) _ x 
x (c ife '- 8 ') Cov(re;0,0') . 



(Dl) 



(D2) 



The subscripts in the exponentials denote real (i, j = 1) 
and imaginary (i, j = 2) parts. The power spectrum P(k) 
is directly related to the covariance of A. In fact, we have 



P(k) = Cov l4 (A; k, k) + (A 4 (fc))(A l (fe)) , 



(D3) 



with summation implied on i and mean over all directions 
of k. For a large set we can perform integrations over 
the whole plane. Thus we find 

P(k) = ^cSv(re;fc) + -l|re (fe)| 2 |pW(fe) - if . (D4) 

Here, as usual, hats indicate Fourier transform. In the 
weak lensing limit the second term of this expression can 
be dropped. Hence, if W is a Gaussian of the form of 
Eq. (3) we find Eq. (47). [Note that no averaging over the 
direction of k is needed in the weak lensing limit if the 
weight function has the form (3).] 
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